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Abstract. Despite the great success of the Huxley sliding filament model proposed half 
P-i a century ago for actin- myosin linkages (cross-bridges), it fails to explain the force- velocity 

6 

5*3 proportions is therefore reconsidered and a new solution to the problem of steady mus- 

O 

cle eccentric and concentric contractions is reported. Instead of numerical modeling the 

>^ 

contractive-force data by appropriate choice of the seemingly arbitrary heterogeneity of at- 
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behavior of stretching skeletal muscles. Huxley's two-state kinetic equation for cross-bridge 



tachment and detachment rates of myosin heads to actin filament cites, Huxley's idea on 
mechanical equilibrium is probed into thermodynamic equilibrium in the whole overlapped 
actin-myosin zone. When the second law of statistical thermodynamics is applied to cross- 
bridge proportions, the weakly bound states appear to be correlated to the strongly bound 



> 

oo 
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states via structural and kinetic intrinsic muscle characteristics. A consequent substantial 
reduction of the number of free parameters in cross-bridge proportions is also due to the 
j>! overall self-consistency (normalization) of attachment-detachment stochastic events. The 

^ explicit force-velocity curve is found to be generic when applied to the reduced tension 

in a single cross bridge, sarcomere, fiber, or muscle as a whole during its active shorten- 
ing or lengthening. This universal curve fits the empirical tension-velocity data on frog 
muscle shortening using only one adjustable parameter, while the Huxley model employed 
four parameters. The established normally distributed cross-bridges, detaching slowly near 
equilibrated states in steady lengthening muscle and quickly in shortening muscle, are in 
qualitative agreement with recent data on the force enhancement following muscle stretching. 
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1. Introduction 

The early studies of muscle fibers under the light microscope revealed cross-striations 
running normal to the fiber axis. It was observed that during either concentric or eccentric 
contractions, the length changes occurred via an increase or decrease in the extent of the 
I-band with the A-band remaining unchanged. Two groups laid the foundations for the cross- 
bridge (CB) theory when they simultaneously suggested that muscle contractions occur due 
to the relative sliding of the thick myosin filaments past the thin actin filaments, mediated 
by the ATP-dependent actin-myosin linkages (H. E. Huxley and Hanson, 1954) working as 
independent force generators (A. F. Huxley and Niedergerke, 1954). 

A. F. Huxley (1957) evaluated muscle tension caused by shortening, in fact, based on the 
idea of the existence of mechanical equilibrium of myosin heads at the regular sites of actin 
filaments. His famous two-state (bound-unbound) sliding filament model was determined on 
the basis of the simplest standard kinetic equation controlled by the velocity-independent 
rates of attachment and detachment of myosin heads. In spite of the great success in illu- 
minating the force generation and power liberation during muscle shortening, the Huxley 
approach generally failed to explain the ascending branch of the phenomenologically estab- 
lished tension- velocity equation (see e.g. Harry et al., 1990 and references therein). 

Exploring the fact that the sliding filament model leaves a free choice of the attachment 
and detachment rate functions, many researchers successfully simulated a range of muscle 
properties in lengthening by fitting the empirical data by linear functions and constants 
suggested by Huxley (1957) for the CB proportions or by bilinear and exponential functions. 
Likewise, considerable efforts have been made to modify the rate functions (Zahalak, 1981; 
Harry et al., 1990; Ma and Zahalak, 1991; Cole et al., 1996) or to find an exact numerical 
solution to Huxley's model (Wu et al., 1997). Very recently, controversies surrounding 
Huxley's approach were brought forth by Mehta and Herzog (2008) in their careful studies 
of force exposed by a single CB during lengthening. 

The theoretical problem of self-consistency in the two-state sliding filament models was 
thoughtfully discussed by Hill and co-workers (Hill et al., 1975). Considering the conditions 
of CB thermodynamic equilibrium besides the minimum of mechanical energy (Huxley, 1957), 
they demonstrated that the original Huxley model has a low efficiency in comparison to its 
modified versions. Moreover, it was noted by Eisenberg et al. (1980) that "there is no in 
vitro evidence for ... the basic (Huxley's) assumption that the cross-bridge detaches slowly". 
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The lacking data were provided by Mehta and Herzog (2008). 

In this study, I develop a statistical approach to the filament sliding mechanism and 
show that only linear functions for the attachment-detachment rates are compatible with 
the concept of thermodynamic equilibrium. A new analytical solution to Huxley's two-state 
kinetic equation is proposed and verified using the available from the literature data on 
tension in steady muscle lengthening and shortening. 



2. Methods 

2.1. Model by Huxley (1957) revisited 

At a fixed muscle contraction velocity V, the number of bound CB states Ny combining 
myosin filament with actin filament of the total number of sites Nqa obeys the common 
"balance" kinetic equation 

d » T , s dNy dNy dx W , T , T . , S „ T 

-^t) = + = f(x)(N 0A -N v )- g ( x )N v . (1) 

Here / and g are attachment and detachment rates of the corresponding unbound and bound 
states located in time t at a distance x estimated from the nearest site x — 0. The steady 
process determined by late times t ^> f~ l , g^ 1 providing dN v (x,oo)/dt = in Eq. (Jl]), 
reduces Eq. (jl]) to 

~^-j-n v (x)= f(x)(l-n v )-g(x)n v , (2) 

i.e., to Huxley's Eq. (4) where the proportion ny(x) = Nv(x,oo)/N a of CBs during 
steady shortening. According to Huxley, the force output F(x) = kx is produced when x 
decreases at a positive velocity of sliding of the actin filament Va and a negative velocity 
of myosin filament Vm, i-e., Va = — Vm = —dx/dt > 0. The contraction velocity per one- 
half sarcomere V/2 determines the contraction velocity V of the muscle as a whole, when 
modeled by V = Va — Vm = 2V \. As can be derived from Huxley's Eq. (6) with the 
preservation in part its notations, the overall generated force 

^(total) ,. sN 0M r L „, , , v dx 



pvan, = Um . F (x)n v (x) — , (3) 

L~>oo L J_ L Zl A 

was evaluated via the force F(x) per one myosin site, as one actin site is carried past it. 
Here s is the sarcomere length, I a is the trial distance between the nearest sites in the actin 
filament evaluated in Huxley's Eq. (15), and Nqm is the number of sites in the thick filament 
in the overlapping zone of length L. 



The solution to Eq. ^ for the shortening regime (hereafter distinguished by index 1) 
was found as a combination of the localized (short-domain) and delocalized (large-domain) 
spatially correlated (bound) states. The corresponding proportions reproduced exactly from 
Huxley's Eqs. (7) and (8) are 
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and Vi = h(fi + gi). 
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In turn, this description of the two CB states follows from the rates postulated by linear 
functions, namely 
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and two constants f'(x) = 0, g'(x) = g[, for — oo < x < 0. The muscle concentric steady 
tension Py reduced to the model isometric tension P found on the basis of Eqs. (j3])-([7]), 
namely 



-,(short) 
V 



V 



exp 



V 
V 



1 + 



VVi 



for V > 0, 



(8) 

10.0 s-\ 



was fitted by the widely cited four model parameters: fi = 43.3 s -1 and gi 
indicating slow detachment of the localized CBs, and f[ = with g[ = 209 for delo- 
calized states. In addition, two more adjustable parameters h « 15 nm and Vi = V^/A, 
where Vm^ is the empirical maximum shortening velocity, were indirectly employed when 
tested by Hill's empirical equation (see Chapter IV in Huxley, 1957). It is noteworthy that 
the nearest-site distance in the actin filament treated as a free parameter was estimated as 
/ w h, i.e. close to the known nearest-molecular distance in the myosin filament Im = 14.5 
nm (Craig and Woodhead, 2006). However, the ratio P_oo/Po = (/i + <7i)/<7i = 5.33 reported 
by Huxley (1957) for the muscle lengthening regime, contrasts to the observed ratios falling 
between 1.8 and 2.0 (e.g. Harry et al., 1990). 



2.2. A new solution to Huxley's steady equation 
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Beyond any specific suggestions, the formal solution to the steady-state Eq. (|2]) 

/ i r x \ dx V 

n v (x) = n (x) + An v (x) = n + (l-n )c v exp f -- / [f(x')+g(x'))dx'\, x = — = =f — , 

(9) 

is valid for any contraction shortening velocity V (= — 2x > 0) and lengthening velocity V 
(= 2a; < 0), leaving an arbitrary choice of the rate functions f(x) and g(x). A differential 
equation of the first order possesses as common only one free constant, denoted by cy, 
whereas 

= mrm (10) 

straightforwardly following from Eq. (J2| taken at V = 0, describes maximal CB pro- 
portions limited by intrinsic rates. In the Huxley model, the constant cy = — noi(l — 
n oi) _1 exp(— V\/V) in Eq. ^ results from his boundary condition ny[h) = providing the 
non- Gaussian proportion Q for CB localized states during muscle shortening. 

Besides the boundary conditions considered below, let us employ the property of peri- 
odicity in the overlapping part of the actin filament of length Nd having N occupied cells. 
Since the Huxley proportion nv(x) in Eq. ^ plays the role of the late-time probability of 
finding one of the two myosin heads attached at a position x between two nearest equivalent 
sites (see also Hill et al. 1975, p. 346), the total force output in a finite overlapped zone is 

i ^ f +Nd dx' f d dx 

F (zone) = NFv= F(x')ny(x / )— ^- — N F{x)n v {x) -J, where x' = xN. (11) 
J-Nd 2rf J_ d 2d 

Here F(x) is the active force per one actin site, substituting that per one myosin cite in Eq. 
(|3|. Such a consideration suggests the statistical equivalence of all the occupied cells in the 
actin filament treated as a one-dimensional crystal of lattice constant d (= 36 nm, e.g. Hill 
et al., 1975). In order to be consistent with Eqs. (11) and ([3]), the normalization conditions 
for the CB distributions (proportions) 

L nv{x) 2% = I nv{x)d i = L Mx)d i = 1 (i2) 

must be taken into consideration. 

Furthermore, extending Huxley's idea on the minimum of mechanical energy at x = over 
the minimum of Gibbs energy (Eisenberg et al., 1980), the CB state with F(0) = is treated 
as the locally equilibrated state, having maximum configurational entropy at x — 0. As the 
consequence of one of the most general principle (second law) of statistical thermodynamics, 



the distribution of bound states ny(x) given in Eq. ([9]) must have Gaussian form centered 
at x = (see e.g. Chapter 12 in Landau and Lifshitz, 1989). One can see from Eq. 
that the thermodynamical principle ensured by the sign requirement xx > can be satisfied 
solely by the linear parameterization of the rate functions, namely 



f{x) = f m —, d( X ) = 9m — 
■Em. *Em. 



for x m = x + > x > or x m 



-X- < x < 0, 
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(14) 



readily follows from the normalization conditions (12), where the standard error function 



erf(y) = (2/y / 7r) J Q y exp(— t 2 )dt, lying between [= erf(0)] and 1 [= erf(oo)], is employed. 



3. Results 

3.1. CB proportions in steady muscle shortening and lengthening 

As seen in Eq. (|9]), a given CB is characterized by the equilibrated velocity-independent 
ground state and the excited non-equilibrated state described, respectively, by uniform pro- 
portion no = fm/(fm + 9m) and non-uniform, heterogenous stochastic proportion Any(x), 



having the meaning of probabilities normalized in Eq. (12). The requirement of signs ( 
xx > 0), while in particular ensuring the self-consistency with the ground state (when 
x — > 0, Any(x) —>■ 0), also constrains possible domains for both CB states, as shown in Eq. 



(13). Indeed, both kinds of domains (0 < x < xq and — xo < x < 0, otherwise no(x) = 0) are 
generally possible for the ground state, whereas only the negative domain, — x_ < x < 0, 
satisfies the requirement of self-consistency of the solution given in Eq. Q. In this way, the 
CB boundary conditions are not postulated as in Eqs. Q-([7]), but result from the minimum 
of Gibbs energy, also giving rise to the CB mechanical constraints, consistent with the si- 
multaneous observations of directions of both the force output and contraction velocity, as 
illustrated in Fig. 1. 

Place Fig. 1 
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The analysis in Fig. 1 specifies domains of the short- domain CB states which, being 
incorporated in the trial Eq. ^ with the help of Eq. (14), yield 



Xq d 2 ex P [~v^~ 

n v (x) = n Q (x) + (1 - n -j) 7= , " T / B v (x). (15) 

a X- y/irv er f ( J_ 



Here, the auxiliary functions Qq(x) = 0(±x) — Q(x =F xq) and Gy(x) = 0(— x) — Q(x + x_) 
are introduced by the standard Heaviside (step) function Q(y), which is one for y > and 
zero for y < 0. 

Muscle tension in steady shortening and lengthening 
The mean force output Fy generated by a single cell of the actin filament is evaluated 



using Eqs. (11) and (15), namely 



_ = F - k — 1 - n -r F = ± — 
2d 2 v d ' v ; 2d 



F y = F + AF y = fc / xn v (x)— = F - fc— (1 - n -^)$(w), F = ±-^> , (16) 



via the CB stiffness k = F(x)/x (Huxley, 1957; Huxley and Simmons, 1971), shown to be a 
velocity- independent intrinsic muscle quantity (e.g. Lombardi and Piazzesi, 1990), where 

= 2<[^ 1 ~?^T^ ; with $(0) = 0, $(1) = 0.846, and $(00) = 1. (17) 

V n erf (^) 



The upper and lower signs in the velocity-independent limiting steady force F (16) cor- 
respond to shortening and lengthening (see Fig. 1). In this way, Huxley's Eq. ^ is 
transformed into a unique equation 



p v Fv . , .,, (d-n x )x^ 

— = — = IT o- m ® (v) , with cr m = 2 > ( 18 ) 

-n -to ^0^0 



for the reduced CB tension and force output in both concentric and eccentric muscle con- 



tractions conducted at positive and negative steady velocities V = vV m (14), respectively. 



The one-parameter fitting analysis of the proposed theory is conducted on the basis of 



Eq. (18) and the available experimental data. In Fig. 2, the muscle shortening is described 
by 

p (short) $ ( ^Y- 

Poi $(A) ' S S max ' 

Knax AV^i, Vml Xml\fml 9ml) i (19) 
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where A is an adjustable parameter. 



Place Fig. 2 



Place Fig. 3 



The steady muscle lengthening is fitted in Fig. 3 by 

p(stret) / V \ 

v - 1 + $ ( — - 1 , -OO < V < 



^02 V ^m2 , 

V m 2 = -X m2 {f m 2 + QnO) < 0. (20) 

using the characteristic velocity V m 2 as a free parameter. Other parameters describing 
two distinct regimes are specified as x_ = x m i, f m = f m i, g m = g ml , for shortening, 
and x_ = x m2 , f m = fms, 9m = 9m2, for lengthening. One can see that V max plays the 
role of the maximum shortening velocity at which p^ hort ^ = o ? and V m 2 is a characteristic 
velocity separating slow and fast lengthening. Also, the limiting tension in the fastest steady 
lengthening is = 2P . 



3.3. CB domains 

The force- velocity fitting analysis alone does not provide details on the CB attachment- 
detachment rates or the domains. Physically, these domains follow from the conditions of 
realization of thermodynamic stability described by a minimum of Gibbs energy (Hill et al., 
1975). Nevertheless, the overall curve conditions of observation can be established here by 



the inequalities cr^ p ^ > 1 > c^ p \ resulting from Eqs. (19) and (20), where the fitting 
parameter <J^ P ^ = $(0.85) _1 = 1.22 is found for muscle shortening and o~^ p \ generally 
lying between 0.8 and 1.0, for lengthening. 

Alternatively, the observation conditions of the predicted branches of the master curve 
can be reformulated in terms of the CB rigor state proportions noi < ^mi^oi (^oi +i m i) _1 



and no2 > a^cfe^ 1 ( x 02 + x^) -1 , obtained with the help of Eq. (18). This finding can 
be improved when the CB geometrical constraints shown in Fig. 1 are taken into account. 
Indeed, since the tail of the myosin molecule is longer than heads, one should expect a 
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geometrical constraint x m 2 > x m i, providing no2 > n>oi- Under the simplified requirements 
of periodicity (x 01 + x ml = d and x 02 = x m2 = d), the CB proportions underlying the 
observation of the master curve are specified in the insets in Figs. 2 and 3. 

4. Discussion 

Huxley's model of the establishment of mechanical equilibrium of myosin heads near 
actin-filament sites is based on the simplest kinetic equation determining a balance be- 
tween unbound and bound actin-myosin states. In a muscle contracting at constant velocity 
V, these two states are described by the proportions 1 — n v (x) and n v (x), satisfying the 
steady-state kinetic equation §2§ at generally arbitrary rates f(x) and g(x). Such a property, 
following evidently from the solution ny (x) found for a general case in Eq. ^ , implies that 
the kinetic equation accounts for the most general features of muscle relaxation, regardless of 
details underlying the attachment-detachment mechanism of myosin heads. Consequently, 
theoretical studies exploring an arbitrary choice of the functional form of the attachment- 
detachment rates, which involved increasing number of numerical parameters, guarantee 
good fit to phenomenological data, but do not shed light on the muscle intrinsic character- 
istics. 

After work by Rayment et al. (1993) on the structural study of force generators in 
contracting muscles, the observations of catalytic domains of myosin being initially weakly 
attached to actin are commonly associated with the weakly bound CB states, and the fol- 
lowing structural changes resulting in tight binding of actin-myosin linkages are associated 
with strongly bound CB states. Since the steady-state equation ^ is a late-time part of 
more general kinetic equation 0, the proportions are also part of the non-steady 

solutions, as demonstrated by Lombardi and Piazzesi (1990) and recently by Walcott and 
Herzog (2008) employing Huxley's Eqs. Q and It seems therefore plausible to associate 
Huxley's short-domain proportion n v ° c \x) and large-domain proportion ny doc \x) with re- 
spectively weak and strong late-time CB states. In this study, the actin-myosin bound state 
is composed of the equilibrated and excited states distributed by Gaussian function dictated 
by the second law of thermodynamics. 

Within the proposed framework of stochastic approach to the attachment-detachment 
events of myosin heads, a common requirement of normalization of the random proportion 
ny (x) specifies the heterogeneity of the CB distribution via the correlated intrinsic structural 
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(x m , d), kinetic (f m , g m ) and dynamic (V m ) muscle characteristics [see e.g. Eq. (14)], that 
decreases the number of free parameters. Moreover, the trend of weakly bound myosin 
heads to achieve maximum structural-domain entropy in the vicinity of actin-filament sites 
(x ~ 0) requires a correlation in signs between the head displacements (x) and velocities (x). 
Consequently, the conceivable CB domains for both bound states schematically shown in 



Fig. 1 are eventually described by Heaviside functions in Eq. (15). A geometrical selection 
of the main components of the force resulting in the power stroke in a direction consistent 
with the vector of contraction velocity are also shown in Fig.l. 

The explicit solution (J9| to Huxley's kinetic equation for CB proportions n (x) and 
Any (a;), distributing respectively strongly and weakly bound states over the actin filament 
cells, results in the velocity-independent (isometric) force Fq and contractive force AFy, 



components of the CB force output Fy (16). In Fig. 2, famous Huxley's comparative 
analysis with Hill's data on muscle concentric tension (Huxley, 1957, p. 287) is revisited. 
The high-velocity wing of the tension curve above V/V max = 0.5 controlled mostly by weak 
CB states is well fitted by both non-Gaussian rt5J) and Gaussian (15) proportions. It is 
not the case of the low velocity region 0.2 < V/V max < 0.5, where a discrepancy between 
Huxley's curve (!§]) and the data indicate a disadvantage of the short-domain (x < h < d) 
weak CBs exerting negative force AFy and eventually reducing the total produced tension. 
In contrast to the postulated retarded detachment (g 1 < fx) discussed in Eqs. (HI) and (]7L 
the Gaussian strong and weak CBs require faster detachment than attachment {g m i > fmi) 
in muscle shortening, as derived from Hill's data and shown in the inset in Fig. 2. The 
regular deviation of Gaussian CBs from the data at very low velocities is associated with 
a simplified modeling of the channel of relaxation of weak states to strong states. Indeed, 
the fit analyses can be improved when the proportion of a new weak-to-strong transient 
Gaussian CB spreads its domain symmetrically within the range —5 < x < 5, extending the 
weak CB state in the vicinity of x ~ 0, as shown by the dotted line in the inset in Fig. 2 
for the case V/V max = 0.1. 



In Fig. 3, the upper branch of the force- velocity curve (18) drawn at a single adjustable 
parameter (V m 2 = —288 nm/s) fits well the empirical data on CB force in muscle stretching. 
Similar to shortening, the fitting analysis could be improved at low stretch velocities when a 
transient bound state is additionally introduced, as independently proposed by Mehta and 
Herzog (2008, Fig. 3). These authors also raised the central question on the existence of 
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Huxley's proportions favoring myosin head attachment events at large distances with an in- 
crease in contraction velocity. One therefore infers that although non-Gaussian large-domain 
CBs ([5]) numerically fit the empirical data (Fig. 2), they do favor neither thermodynamic 
equilibrium in the overlapped zone nor high cycle efficiency (Hill et al., 1975). 

The microscopic structures of the short-domain bound states are provided above via ob- 



servation conditions of the generic curve (18) equally applied to the reduced tension in a 
single CB, sarcomere, fiber, or muscle as a whole during its steady shortening or length- 
ening. It is also demonstrated (inset in Fig. 3) how the two-state muscle cycle duty ratio 
ft derived from real experiments can be helpful in a characterization of the Gaussian CB 
rate ratio a = f m /g m [= (1 — 0)/j3\ and strong bound state proportion n (= 1 — (3). To 
summarize a comparative analysis of structural and kinetic characteristics of CBs, one can 
see that steady muscle eccentric and concentric contractions are well distinguished via the 
attachment-detachment rate rations, with a stret > 1 > a short , the strongly-bound occupa- 
tion numbers, with n^ tret ^ > n^ hort \ supported by the directly observable cycle duty ratios 
Pstret < Pshort (Mehta and Herzog, 2008). Within this context, the working hypothesis by 
Mehta and Herzog (2008) "that a stretched cross-bridge might remain attached longer than 
a cross-bridge that had been shortened while attached" combined with the main finding by 
Lombardi and Piazzesi (1990) that "reattachment (in steady lengthening is)... faster than 
attachment in the isometric condition or during shortening ... in the same domain of x" 
results in the predictions a stre t > 1, n^ tret ' > > 1/2, and f3 stret < 1/2, which are generally 
consistent with the CB parameters derived in the insets in Figs. 2 and 3. 

To conclude, the provided statistical thermodynamic analysis of the attachment- 
detachment CB process, modifying Huxley's mechanical sliding filament model, can also 
be figured out as an two-headed steady walking of synchronous myosin molecules over peri- 
odical sites of actin filaments with multiple 36-nm steps, as directly observed by Sakamoto 
et al. (2008). The fluctuating steps are statistically scattered by the normal distribution, 
having the zero mean and variance linear with muscle contraction velocity. The proposed 
steady contraction dynamics is universally observable through the two branches of the force- 
velocity curve generic for steady shortening and lengthening of a muscle as a whole or its 
counterparts. The microscopic structural muscle characteristics appear to be strongly cor- 
related to kinetic and dynamic characteristics distinguished by the force output directions 
generated in distinct muscle regimes. At a macroscopic level, similar kind of correlations 
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driven by maximum generated force were revealed via the primary muscle functions well 
distinguished though the muscle structure adapted to efficient eccentric, isometric, or con- 
centric contractions (Kokshenev, 2008). 
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Figure Legends 



Shortening, V - - 2dx/dt >0 Lengthening, V = 2dx/dt < 




Fig. 1 

Figure 1. Mechanical scheme of the force generation by combining myosin heads with 
periodic actin filament. Each of the two heads of the effective CB may be attached to actin 
filament either in equilibrated ground state (shown by the open circle) with the uniform 
probability n , within the domains x < ±x , or in the non-equilibrium, excited state (closed 
circle) with the probability Any(x), within the domains x < —xi, —x<i- The arrows indicate 
the directions of the sliding velocity of the actin filament Va and the myosin filament Vm- 
During concentric muscle contraction with a positive velocity V, the velocity-independent 
portion of the generated force F is also positive, whereas the ATP hydrolysis results in the 
negative portion of the contractive force AF V . During eccentric contractions commonly 
associated with the negative direction of velocity V, both the forces are also negative. 
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Figure 2. Analysis of the theoretically predicted tension- velocity curve using available 
data on the reduced tension during steady muscle shortening. The points and dashed- 
point curve are the famous data by Hill (1938) for isolated frog muscles modeled by Huxley 
(1957), drawn respectively by the phenomenological equation p( cxp ) / P = a(l—V/V mgx ) / (a+ 
V/Vmax), with a = 0.25 and Eq. (8), fitted by Huxley's parameters listed above. The solid 
line is Eq. (19) taken at A = 0.85. Inset: The attachment-detachment rates and CB 
proportions predicted in Eq. (15) within the CB domains at distinct shortening velocities 
reduced to the maximum velocity. The CB structure discussed in the Results is exemplified 
by the model parameters x^ od ^ = 2d/3 and x^° d ^ = —d/3, as well as by n^ od ^ = 0.47, 
providing the rate ratio a\ = f m \/g m \ = 0.88. The dotted line shows a proportion for the 
modeled transient CB state schematically drawn for V/V max = 0.1. 
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Figure 3. Steady force induced by one cross-bridge versus the stretching velocity. The 
open circles are the mean datapoints of the forces (re-scaled by |F | = 1.95 pN ) measured by 
Lombardi and Piazzesi (1990, Fig. 7) in frog muscle fibers at stretch velocities lying between 
75 and 1030 nm/s and scaled here by d = 36 nm. The closed square indicates the force 
per one CB reported by Mehta and Herzog (2008) for unspecified velocities. The theoretical 
curve is drawn based on Eq. (20) with V^° d ^ = —8 dj s. Inset: The attachment-detachment 
rates within the CB domains and CB proportions predicted in Eq. (15) at three distinct 
velocities reduced to the found V^ oA \ They are exemplified by model parameters x^°^ = 
x m2^ = consistent with the observation conditions discussed in the Results, as well as 
by = 1 — /?2 CXP ' ) = 0.93, where the stretch cycle duty ratio (3^ xp ^ = 7.35% (the time of 

attachment J" 1 related to total CB cycling time f^+g^, i.e., (3 = 1 — n ) studied by Mehta 
and Herzog (2008) is employed. Moreover, the relation |V^™ od ' ) | = fm^^m^ /[I — fii^} 
derived from Eq. (20) provides a crude model estimate for the CB domain ~ 1-2 d, if 

their characteristic attachment time [/a^] -1 = 0.167 s is also employed. The rates (shown 
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by dashed lines) are determined by the ratio / 9m2 = 13, corresponding to the model 
estimate [g^?^] -1 = 2.2 s. The dotted line shows a proportion for the assumed transient 
CB state schematically drawn for V/V^ od ^ = 0.1. 
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